Accepted for publication in Solar Physics. 

The authoritative version is DOI: 10.1007/sll207-012-0149-8 



Resolving the Azimuthal Ambiguity in Vector 
Magnetogram Data with the Divergence- Free 
Condition: the Effects of Noise and Limited Spatial 
Resolution 

A.D. Crouch 

Received: 6 May 2012 / Accepted: 24 September 2012 

Abstract We investigate how the azimuthal ambiguity in solar vector magne- 
togram data can be resolved by using the divergence-free property of magnetic 
fields. In a previous article, by Crouch, Barnes, and Leka {Solar Phys. 260, 271, 
2009), error- free synthetic data were used to test several methods that each make 
a different assumption about how the divergence-free property can be used to 
resolve the ambiguity. In this paper this testing is continued with an examination 
of the effects of Poisson photon noise and limited instrumental spatial resolution. 
We find that all ctirrently available methods based on the divergence-free prop- 
erty can produce undesirable results when photon noise or unresolved structure 
are present in the data. We perform a series of experiments aimed at improving 
the performance of the global minimisation method, which is the most promising 
of the methods. We present a two-step approach that produces reasonable results 
in tests using synthetic data. The first step of this approach involves the global 
minimisation of a combination of the absolute value of the approximation for 
the divergence and a smoothness constraint, which is designed to minimise the 
difference between the magnetic field in neighbouring pixels. In the second step, 
pixels with measurements known to be strongly affected by photon noise are 
revisited with a smoothing algorithm that also seeks to minimise the difference 
between the magnetic field in neighbouring pixels. 

Keywords: Sun: magnetic field 
1. Introduction 

Knowledge of solar magnetic structures and solar activity in general is greatly 
enhanced by measurements of the vector magnetic field in the solar atmosphere. 
However, there is an ambiguity of 180° in the direction of the inferred component 
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of the magnetic field perpendicular to the line-of-sight, when it is derived from 
the linear polarisation of magnetically sensitive spectral lines (Harvey, 1969). 
Many algorithms are available to resolve this ambiguity in single-height vector 
magnctogram data, and while some methods perform better than others (for 
overviews and comparisons sec Mctcalf et al, 2006; Leka et al., 2009), it is appar- 
ent that each method involves at least one major assumption or approximation 
that may not be valid for solar magnetic fields in general. 

A class of methods that can potentially lessen the need for additional assump- 
tions resolve the azimuthal ambiguity by using the divergence-free property of 
magnetic fields and multiple- height vector magnctogram data {e.g., Wu and 
Ai, 1990; Cuperman, Li, and Scmel, 1993; Li, Cuperman, and Scmel, 1993; 
Boulmezaoud and Amari, 1999; Li, Amari, and Fan, 2007; Crouch and Barnes, 
2008; Crouch, Barnes, and Leka, 2009). In this context, multiple-height vector 
magnctogram data include information regarding the; variation of the magnetic 
field in the horizontal heliographic directions and along the line-of-sight direction 
(the latter of which can be inferred from observations, see, e.g., Ruiz Cobo and 
del Toro Iniesta, 1992; Collados et al, 1994; Metcalf et al., 1995; del Toro Iniesta 
and Ruiz Cobo, 1996; Liu et al, 1996; Westendorp Plaza et al, 1998, 2001; 
Socas-Navarro, Trujillo Bueno, and Ruiz Cobo, 2000; Elbe et ai, 2002; Leka and 
Metcalf, 2003; Socas-Navarro, 2005, 2007). It is important to emphasise that 
the horizontal heliographic directions are not orthogonal to the line-of-sight 
direction, except at disk centre. The problem of computing the divergence in 
terms of horizontal heliographic and line-of-sight derivatives, for any position on 
the solar disk, was addressed by Crouch and Barnes (2008) and Crouch, Barnes, 
and Leka (2009, henceforth Paper I). 

In Paper I, using error-free synthetic data with measurements available at two 
discrete heights, three methods were tested that each make different assumptions 
about how to use the divergence-free property to resolve the ambiguity. The 
synthetic data employed in Paper I did not account for noise that is typical 
of solar observational data. In this article we use more realistic synthetic data 
to test how the currently available methods respond to the effects of noise. 
Following the procedure of Leka et al. (2009) two types of noise are included 
in the synthetic data: noise to simulate Poisson photon noise in the observed 
polarization spectra, and a spatial binning to simulate the effects of limited 
instrumental spatial resohition. 

The outline of this article is as follows. In Section 2 we re-iterate the derivation 
of the divergence-free condition for any position on the solar disk in terms of 
observable quantities, as described in Crouch and Barnes (2008) and Paper I. In 
Section 3 we describe the synthetic data and metrics used to test the performance 
of the ambiguity-resolution algorithms. In Section 4 we test the various methods 
based on the divergence-free property that were examined in Paper I. We find 
that all currently available methods can produce undesirable results when the 
data contain either photon noise or unresolved structure. Nevertheless, we find 
the same general trend that was found in Paper I, in that the global minimi- 
sation method is more robust than both Wu and Ai's (1990) criterion and the 
sequential minimisation method (Boulmezaoud and Amari, 1999). Consequently, 
in Section 5 we continue the development of the global minimisation method by 
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implementing and testing several techniques to minimise the effects of noise on 
the disambiguation results. In Section 6 we draw conclusions. 

2. The Divergence-Free Condition 

In this section we re-iterate how the divergence of the field can be expressed in 
terms of observable quantities for any position on the solar disk, as described by 
Crouch and Barnes (2008) and Paper I. Within a limited field of view we assume 
that a layer of constant optical depth can be represented by the heliographic 
plane, which is tangent to the solar surface at the position given by the central 
meridian angle, the latitude and the P- and _Bo-angles. On the heliographic 
plane, the image and heliographic components of the magnetic field are related 

by 

= aiiB^. + ai2By + aisB], , 

B^ = a2iBl + a22Bl + a23Bi, (1) 

= (Xai-Bi + ^32-8^ + 033^^ , 

where the coefScients aij are given by Equation (1) of Gary and Hagyard (1990), 
the superscipts h and i label heliographic and image components, respectively, 
Bl, = is the line-of-sight component of the field, and the image components 
of the field in the directions perpendicular to the line-of-sight direction are 

bI = B± cos^ , and B'y = B± sin^ , (2) 

where B± is the magnitude of the transverse component of the magnetic field 
(perpendicular to the line-of-sight) and ^ is the azimuthal angle measured coun- 
terclockwise from the a;'-axis to the transverse component of the field (where 
x' and are coordinates on the plane perpendicular to the line-of-sight). In 
solar vector magnetogram data, the azimuthal ambiguity occurs because ^ can 
only be inferred within the range < ^ < 180°, using the linear polarisation of 
magnetically sensitive spectral lines. 

When written in heliographic coordinates, the calculation of the divergence 
requires the derivative of the magnetic field in the direction perpendicular to the 
heliographic plane, d/dz^, which cannot be directly measured away from disk 
centre. However, the derivative of the magnetic field in the line-of-sight direction, 
d/dz\ can be measured for any position on the solar disk in principle {e.g., Ruiz 
Cobo and del Toro Iniesta, 1992; Collados et al, 1994; Metcalf et ai, 1995; del 
Toro Iniesta and Ruiz Cobo, 1996; Liu et al., 1996; Westendorp Plaza et al., 
1998, 2001; Socas-Navarro, Trujillo Bueno, and Ruiz Cobo, 2000; Elbe et al., 
2002; Leka and Metcalf, 2003; Socas-Navarro, 2005, 2007). For more detailed 
discussion of the various methods for inferring the line-of-sight variation of the 
magnetic field see Crouch and Barnes (2008) and Paper I. For all of the available 
methods the magnetic field is inferred at discrete heights (in the line-of-sight 
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direction) and, therefore, the hne-of-sight derivatives of the magnetic field must 
be approximated. The derivatives d/dz^ and d/dz^ are related by 

df df df df 

5^=«13a^+«23^+a33^, (3) 

where / is any differentiable function, and a;^ and -y^ are coordinates on the 
hcliographic plane. The derivatives d/dx^ and d/dy^ can be approximated using 
discrete measurements over the heliographic plane. We use line-of-sight distance 
z\ but it should be noted that the magnetic field may be inferred as a function 
of optical depth t; in such cases, to determine the correspondence between line- 
of-sight distance and optical depth, an atmosphere model would need to be 
employed {e.g., Maltby et al., 1986; Vernazza, Avrett, and Loeser, 1981; Collados 
et al, 1994; Socas-Navarro, 2007). 

Using Equations (1) and (3), the divergence of the magnetic field can be 
expressed in terms of observable quantities {i.e., derivatives of the image com- 
ponents of the magnetic field with respect to x^, and z^), 

a33V-B = £>„ + a33^, (4) 

where 



Da = a3i^-Y+a32^-f + (aiia33- 013031) + (0112033-013032) 



dz' ' dz' ' feh ' ^^^^-^^ ^^^"•^^^ dx^ 

+ (021033 - 023031) + (022033 - 023032) , (5) 

which demonstrates that the calculation of the divergence generally involves 
the line-of-sight derivatives of all three image components of the magnetic field, 
except at disk centre where Uij — dij and only dB^/dz^ is required (see Crouch 
and Barnes (2008) for further discussion). 

It is important to emphasise that the line-of-sight derivatives of the transverse 
image components of the field, dB^^/dz^ and dBy/dz\ depend on the ambigu- 
ity resolution at all heights used to approximate these derivatives. Therefore, 
the ambiguity resolution at one height can afi^ect the ambiguity resolution at 
other heights, when the ambiguity is resolved using an algorithm based on the 
divergence for data positioned away from disk centre. On the other hand, for 
data positioned at disk centre and algorithms based on the divergence, each 
height can be disambiguated independently if the approximation for horizontal 
heliographic derivatives involves only measurements at one height. 



3. Synthetic Data and Performance Metrics 

Leka et al. (2009) compared several techniques for disambiguating single-height 
magnetogram data using synthetic data for which information regarding the 
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variation of the magnetic field out of the heliographic plane was not available. 
They constructed two families of synthetic data that included two different 
types of noise that are expected in solar observational data: noise to simulate 
Poisson photon noise in the observed polarization spectra, and a spatial binning 
to simulate the effects of limited instrumental spatial resolution in the directions 
perpendicular to the line-of-sight. Methods based on the divergence require mag- 
netic field data from at least two heights in order to approximate the line-of-sight 
derivatives of the field. To test the performance of the various algorithms based 
on the divergence we use the same magnetic field configurations and the same 
types of noise as were considered by Leka et al. (2009) with magnetic field data 
available at two heights. The pixels at each height are aligned in the line-of- 
sight direction and separated by a line-of-sight distance A^'. All of the synthetic 
datasets used in this article were provided by K.D. Leka and Graham Barnes 
(2012, private communication). In the following paragraphs we present a brief 
summary of the procedures used to generate the synthetic data. 

To simulate the influence of noise on the observed spectra the two families 
of synthetic data are generated using the same general approach, with the 
following steps (complete details are provided in Leka et al., 2009): (1) At 
each pixel at each measurement height an analytic magnetic field is computed 
(described below). (2) Unno-Rachkovsky-Stokes profiles are generated using a 
Milne-Eddington atmosphere, with parameters consistent with the A630.25 nm 
Fe I line, and subsequently convolved with an instrumental response function. 
(3) The resulting Stokes profiles are modified according to the specifics of the 
type of noise being tested. (4) The modified "noise-added" Stokes profiles are 
inverted using the same Milne-Eddington algorithm used in step (2). 

For the case with photon noise, the magnetic field is constructed from a 
collection of point sources located on a plane below the heliographic plane. 
The magnetic field is linear force-free and the contribution to the field from 
each source is determined by a Green's function {e.g.. Chiu and Hilton, 1977). 
This magnetogram is located away from disk centre, centred at (9°S, 36°E). To 
mimic the effect of noise in photon-counting instruments, randomly generated 
Poisson-distributed noise is incorporated into the Stokes profiles at step (3). 
Three different levels of noise are added: one no-noise case (skipping steps (2) 
through (4)) and two noise-added cases to simulate noise levels expected in 
commonly used instruments. The separation of the two measurement hciights in 
the line-of-sight direction, A^', is equal to the pixel size. The addition of this 
photon noise causes the inferred magnetic field to deviate from the original, 
no-noise, analytic magnetic field. For the two noise-added cases, the ambiguity- 
resolved "answer" at a given location is taken as the realisation of the transverse 
component of the noise-added magnetic field that is closest to the transverse 
component of the original, no-noise, analytic magnetic field. Figure 1 highlights 
the effects of this photon noise by showing maps of the angular difference between 
the ambiguity-resolved, transverse components of the magnetic field for the no- 
noisc and noise- added cases. Further discussion of this test case can be found in 
Section 3.1 of Leka et al. (2009). 

For the case of limited spatial resolution, a current-free (potential) magnetic 
field is constructed that has fine-scale structure with properties similar to both 
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Figure 1. (a) Angular difference between the transverse components of the magnetic field for 
the no-noise and low-noise cases. For the low-noise case, the ambiguity is resolved by taking the 
realisation of the transverse component of the magnetic field that is closest to the transverse 
component of the no-noise magnetic field. The grey-scale image is designed to saturate at 45° . 
Positive / negative vertical magnetic flux for the no- noise case is indicated by red /blue contours 
at 50, 700, 1500 and 2800 G, and the corresponding magnetic neutral line is indicated by the 
green contour. This plot is similar to that in the right-hand panel of Figure 3 of Leka et al. 
(2009). (b) Same as (a) except for the high-noise case. 

penumbrae and plage. The model field is computed on a grid with a pixel size 
of 0.03". To simulate the effects of imperfect instrumental spatial resolution 
the Stokes profiles at step (3) are spatially averaged by factors of 5, 10, and 
30, resulting in effective pixel sizes of 0.15", 0.3", and 0.9", respectively. Some 
of the fine-scale features present in the original magnetic field, that are fully 
spatially resolved on the 0.03" grid, are unresolved at lower spatial resolution. 
The separation of the two measurement heights in the line-of-sight direction, A2:', 
is 0.18" for all pixel sizes. In this case, the ambiguity- resolved "answer" at a given 
location is taken as the realisation of the transverse component of the magnetic 
field that is closest to transverse component of the directly binned (spatially 
averaged), original, analytic magnetic field. For some additional discussion of 
this test case see Section 3.2 of Leka et al. (2009), and Leka and Barnes (2012), 
Georgouhs (2012) and Leka et al. (2012). 

For each synthetic dataset, the magnetic field is sampled at (exact) discrete 
spatial locations in the line-of-sight direction. Thus, the effects of limited spatial 
resolution in the line-of-sight direction are not included in these synthetic data. 
The separation of the pixels in the line-of-sight direction, Az', is assumed to 
be constant over the field of view, although this may not be true for solar 
observations depending on how the solar atmosphere varies in the directions 
perpendicular to the line-of-sight. The synthetic datasets include the effects of 
curvature; however, the layer where the magnetic field is measured is assumed to 
be at constant geometrical height, which may also not be accurate for solar ob- 
servations. Clearly, the synthetic datasets employed in this article are designed to 
isolate some, but not all, sources of noise and/or uncertainty that are expected in 
solar observations. The sources of noise that are not modelled in these synthetic 
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datasets may be potentially significant for disambiguation methods that rely on 
line-of-sight derivatives of the magnetic field and, therefore, further investigation 
is required to quantify these effects. 

To quantify the performance of the various algorithms we use the following 
metrics: the fraction of pixels correctly disambiguated A^area, and the fraction of 
transverse field above a specified threshold T correctly disambiguated A4b±>Tj 
for which we present results for two thresholds T: 100 G (gauss) and 500 G 
(Leka et al, 2009). 

4. The Currently Available Algorithms 

In this section we test the three disambiguation methods that were examined 
in Paper I, with each algorithm implemented exactly as described therein. The 
algorithms employ different assumptions and implementation details regarding 
how to use the divergence-free condition to resolve the azimuthal ambiguity. 

4.1. The Wu and Ai (1990) Criterion 

Wu and Ai (1990) outlined an approach that was later tested on synthetic data 
by Cuperman, Li, and Semel (1993), Li, Cuperman, and Semel (1993), Li, Amari, 
and Fan (2007) and Paper L Following Paper I, for any position on the solar disk 
(except the limb), the divergence- free condition can be written as an inequality, 
i.e., 



In principle, Equation (6) can be used to resolve the ambiguity by choosing the 
realisation of the azimuthal angle such that Da satisfies the inequality, given only 
the sign of dB\/dz^, not its magnitude. A couple of limitations of this approach 
were demonstrated in Paper I; these include (1) the results are sensitive to the 
smoothness of the initial configuration of azimuthal angles and the order in which 
pixels are visited; and (2) the underlying assumption made by this method can 
be incorrect over a significant fraction of the field of view when the derivatives are 
approximated by discrete measurements (in that the sign of Da can be incorrect 
for the correct choice of azimuthal angle, or the sign of Da can be the same for 
both the correct and the incorrect choices of azimuthal angle). 

We examine the case where the initial configuration of azimuthal angles is 
such that By > 0. Results for the noise-added test cases are provided in Figure 3 
and Table 1, and those for the tests of limited spatial resolution are provided in 
Figure 4 and Table 2. In the interest of saving space we only show results for 
the lower height of each magnetogram; however, it is worth noting that we find 
the same trend found in Paper I, where the results retrieved by the Wu and Ai 
criterion at the upper height tend to be worse than those at the lower height 
for data positioned away from disk centre (this is because the upper height is 
disambiguated after the lower height). 
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For the multipole field configuration without noise, the results for the Wu and 
Ai criterion (see Figure 3 and Table 1) are similar to those; found in Paper I. 
We find that the assumption made by this method is violated over a significant 
fraction of the field of view (results not shown, but see Figures 4 and 5 of Paper I 
for demonstrations of this issue). Since the underlying assumption is not tested 
during the disambiguation procedure this method appears to do quite well for 
the case without photon noise. However, as the level of photon noise increases 
the results produced by this method deteriorate significantly. 

Regarding the tests of limited spatial resolution (Figure 4 and Table 2), for 
the higher resolution cases the solutions retrieved by the Wu and Ai criterion 
appear to be reasonable; however, we again find that the assumption made by 
this method is violated over a significant fraction of the field of view (for all pixel 
sizes). For the cases with pixel sizes 0.15" and 0.3" the quality of the solutions 
decreases gradually as the spatial resolution degrades. On the other hand, the 
solution retrieved for the lowest resolution case (pixel size 0.9") is noticeably 
worse than the other cases. 

4.2. The Sequential Minimisation Method 

The sequential minimisation method resolves the azimuthal ambiguity by assum- 
ing that the approximation for the magnitude of the divergence at each pixel 
is minimised for the set of pixels used in the approximation. This assumption 
is applied sequentially at each pixel in the field of view. This method was first 
proposed by Boulmezaoud and Amari (1999) and was examined in Paper I. For 
the test cases examined in Paper I it was shown that the sequential minimisation 
method is more robust than the Wu and Ai criterion but it can produce errors 
when the underlying assumption is violated. The underlying assumption is that 
the correct choice of azimuthal angle will be retrieved for the set of pixels used to 
approximate the divergence, provided that any previously resolved pixels have 
the correct choice of azimuthal angle. Moreover, in Paper I it was demonstrated 
that the sequential minimisation method can propagate erroneous solutions into 
regions that do not violate the underlying assumption. This method is not sen- 
sitive to the initial configuration of azimuthal angles although it is sensitive to 
the order in which pixels are examined. 

Results for the sequential minimisation method are provided in Tables 1 and 
2 and Figures 3 and 4. For the multipole field configuration without noise the 
results are worse than the similar case examined in Paper I. This is primarily 
due to the propagation of erroneous solutions. As the degree of photon noise 
increases many more pixels violate the underlying assumption (not shown but 
see Figures 6(b) and 7(c) and 7(d) of Paper I for a demonstration of the issue). 
Consequently, the solutions retrieved by this method get progressively worse as 
the level of noise increases. As in Paper I, for data positioned away from disk 
centre we find that the results retrieved by this method at the upper height (not 
shown) are very similar to those found at the lower height. For the tests of limited 
spatial resolution (Figure 4 and Table 2) we find again that the underlying 
assumption made by this method is violated for a small fraction of pixels (for all 
pixel sizes), and propagation of erroneous solutions accounts for a large fraction 
of the incorrect results retrieved by this method. 
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4.3. The Global Minimisation Method 

The global minimisation method resolves the azimuthal ambiguity by assuming 
that the correct configuration of azimuthal angles corresponds to the minimum 
of 

^ = EE + |(V-B),,,.,+i|) (7) 

where (V-B)jjjt is the approximation for the divergence at pixel {i,i,k). The 
indices i, j and k label pixels in the x'-, y'- and ^'-directions, respectively, and 
Tlx and Uy are the number of pixels in the a;'- and t/'-directions, respectively. 
Simulated annealing {e.g., Metropolis et al, 1953; Kirkpatrick, Gelatt, and Vec- 
chi, 1983; Press et al., 1992) is used to search for the configuration of azimuthal 
angles that corresponds to the global minimum of E. For the error-free synthetic 
data used in Paper I it was demonstrated that the global minimisation method 
generally produced better results than both the Wu and Ai criterion and the 
sequential minimisation method. However, it was shown in Paper I that the 
assumption made by the global minimisation method (that the global minimum 
of E does indeed correspond to the correct configuration of azimuthal angles) 
can be violated when the divergence is approximated from discrete measure- 
ments, depending on how the magnetic field varies and how it is sampled by the 
observations. 

The simulated annealing algorithm is implemented as described in Paper I 

with the following modifications. Here, we use a "temperature" that can be 
position-dependent. By using a temperature that is local this modification im- 
proves the efficiency of the simulated annealing algorithm in comparison to one 
with a uniform temperature. Moreover, an algorithm with a position-dependent 
temperature tends to retrieve solutions with lower values for E (when noise is 
present) than an algorithm with a uniform temperature (all other things equal). 
Starting with the initial configuration of azimuthal angles we examine lOOOn^Uy 
reconfigurations of randomly chosen pixels at each height that is disambiguated 
(which are all accepted) and take two times the maximum value of \AE\ in a 
small neighbourhood as the initial temperature at each pixel, where AE is the 
change in E (Equation (7)) caused by changing the azimuthal angle. 

Because simulated annealing is a stochastic algorithm, we disambiguate each 
synthetic magnetogram 20 times with different sequences of random mimbers for 
each disambiguation. We select the solution that produces the lowest value of 
E for further examination; results are provided in Figures 3 and 4 and Tables 1 
and 2. 

For the multipole field positioned away from disk centre without photon noise, 
the global minimisation method retrieves the correct solution at every pixel over 
both heights, consistent with the findings of Paper I (for cases with comparable 
values for Az'). As the level of photon noise increases the quality of the results 
produced by the global minimisation method decreases substantially (Figure 3 
and Table 1). As found in Paper I, areas with the incorrect ambiguity resolution 
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tend to occur at similar locations at both heights for data positioned away from 
disk centre. Photon noise tends to have the greatest influence on the results 
retrieved by the global minimisation method in regions where the component 
of the magnetic fleld transverse to the line-of-sight is weak (Leka et ai, 2009): 
Note the trend in A1_b^>ioo g for the different noise levels compared to that in 
-^Bj_>500 G and see Figure 1. For both noise-added cases, E for the solution re- 
trieved by the global minimisation algorithm is less than that for the answer: For 
the low-noise case the ratio of E for the retrieved solution to that for the answer 
is 0.94, whereas for the high-noise case the ratio is 0.84. This demonstrates that 
the assumption made by the global minimisation method is violated for these 
test cases. Moreover, the results produced by the global minimisation method, as 
implemented here, are relatively poor when compared with the results produced 
by the better performing single- height algorithms tested in Leka et al. (2009). 




Figure 2. An estimate of tlie degree of unresolved structure in each 0.9" pixel as a function 
of and j/'^. For each 0.9" pixel we examine the underlying original model magnetic field 
projected onto the image plane ( S^, By ), as sampled on the 0.03" grid. Shown is the cosine 
of the largest angular separation between any two vectors ( S^, By ) that contribute to each 
0.9" pixel. Light/dark shaded pixels correspond to areas with a large/small range of angles; 
the range shown is [—1, 1]. Pixels with large values of this quantity are a reasonable predictor 
of the locations where the global minimisation method fails to retrieve the correct azimuthal 
angle. The magnetic neutral line is indicated by the green contour. 

For the tests of limited spatial resolution (Figure 4 and Table 2) we again 
find that the assumption made by the global minimisation method is violated, 
in that E for the retrieved solutions is less than that for the answers. The ratios 
of E for the retrieved solutions to those for the answers are 0.89, 0.94 and 0.98 
for pixel sizes of 0.15", 0.3" and 0.9", respectively. Despite this, for the higher 
spatial resolution cases the solutions retrieved by the global minimisation method 
are quite reasonable for most of the field of view. The exception to this is the 
plage region (centre, right of the field of view), which evidently poses a challenge 
for this method even at the higher spatial resolution. The plage region consists 
of the collection of randomly distributed small-scale azimuth centres with a 
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significant horizontal magnetic field. The strength of the vertical component of 
the magnetic field of each small-scale structure is proportional to cxp(— r^/a^), 
where a is 0.156" and r is the distance from the azimuth centre (Leka et al., 
2009). Thus, the small-scale structures in the plage region are spatially resolved 
on the 0.03" grid, but are spatially unresolved at lower spatial resolution. 

For the lowest resolution case (with pixel size 0.9") incorrect solutions are 
retrieved by the global minimisation method in and around regions with sig- 
nificant unresolved structure (sec Figure 2). These regions include the plage 
region, the area between the two flux concentrations in the bottom, centre of 
the field of view, and around the fiux concentration at the top, right of the field 
of view (where there is a ring of small-scale azimuth centres that are spatially 
unresolved at lower spatial resolution, see e.g., Figures 5 and 7 of Leka et al. 
(2009)). In these locations disambiguation methods are not expected to retrieve 
good solutions as the magnetic field vectors underlying each 0.9" pixel have a 
large range of angular differences (see Figure 2) and this information is lost at 
low spatial resolution. 
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Figure 3. Disambiguation results at the lower height for the various methods applied to the 
multipole field with various degrees of noise. Areas that are correctly/incorrectly resolved are 
black/white. Positive/negative vertical magnetic flux for the noise-free case is indicated by 
red/blue contours at 50, 700, 1500 and 2800 G, and the magnetic neutral line is indicated 
by the green contour. (Left to right) No-noise, low-noise and high-noise cases. (Top row) the 
Wu and Ai (1990) criterion with initial configuration By > 0. (Middle row) The sequential 
minimisation method. (Bottom row) The global minimisation method (Equation (7)). 



5. Further Development of the Global Minimisation Method 

In the previous section we found that aU currently available disambiguation 
methods that are based on the divergence are sensitive to both photon noise 
and limited spatial resolution, if implemented as described in Paper I. Generally 
speaking, we found the same trend that was found in Paper I, in that the global 
minimisation method is more robust than both the Wu and Ai criterion and the 
sequential minimisation method, in the presence of both types of noise. However, 
the performance of the global minimisation method is significantly compromised 
when photon noise or unresolved structure is present in the magnetogram data. 



ms.tex; 11 October 2012; 0:23; p. 13 



14 



A.D. Crouch 




Figure 4. Disambiguation results for the tests of spatial resolution. Areas that are cor- 
rectly/incorrectly resolved are black/white. The green contour is the magnetic neutral line. 
(Left to right) 0.15", 0.3" and 0.9" cases. (Top row) The Wu and Ai (1990) criterion in 
which the initial configuration has By > 0. (Middle row) The sequential minimisation method. 
(Bottom row) The global minimisation method (Equation (7)). 



In this section we present results from a series of experiments aimed at improving 
the performance of the global minimisation method. 

5.1. Approximation of Derivatives for the Divergence 

In this section we discuss the approximations used to compute derivatives for the 
various terms in the divergence. For all three algorithms examined in Section 4, 
derivatives with respect to and at {i,j,k) are approximated with three- 
point finite differences using measurements at (i, j, k), + k) and (i, j' + l, fc), 
except at the boundaries i — Ux and j — Uy where the approximation is changed 
to use only pixels within the field of view. Derivatives with respect to are 
approximated with two-point finite differences using measurements from the 
two heights. Although they are simple, these finite-difference approximations 
are evidently sufficient for the global minimisation method to retrieve reasonable 
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results for synthetic data without noise, depending on the nature of the field and 
the resolution of the observations. On the other hand, this is generally not true 
in regions where the influence of noise is significant. 

In regions with significant unresolved structure (see Figure 2), the measure- 
ment at a given pixel is a complicated average of the magnetic field underlying 
the area covered by the pixel. Consequently, for the tests of limited spatial 
resolution the value of the magnetic field at a given point may not be known 
exactly and variations of the magnetic field in the horizontal directions may not 
be resolved. Under these conditions finite differences and finite elements are not 
expected to produce accurate approximations for derivatives. 

When the influence of photon noise is significant, the measurement of the 
magnetic field may deviate substantially from the underlying, no-noise value 
(see Figure 1). Finite differences are therefore not expected to produce accurate 
approximations for derivatives when the influence of photon noise is significant. 
In Figure 5 we demonstrate how different approximations for the derivative 
dBl./dx^ respond to the effects of photon noise. In Figures 5(a) and (b) we 
show maps of the difference between the approximation for dBl./dx^ for the no- 
noise and low-noise ambiguity-resolved answers. It is evident from these maps 
that areas with large errors in the approximation for this derivative tend to 
correspond to areas where the global minimisation method retrieves poor results 
for the ambiguity resolution in the noise-added test cases (for example, compare 
Figure 5(a) with the lower panels of Figure 3). 

Several additional points are worth noting regarding Figure 5: (1) Finite- 
difference and finite-element approximations for derivatives can increase (or 
amplify) the effect of photon noise. (2) We find similar results to those shown 
in Figure 5 for the approximation for the derivative d/dy^ and for the other 
components of the magnetic field (not shown). (3) We find that the discrepancy 
between the approximations for the derivatives in the noise-added and noise- 
free cases generally increases as the level of photon noise increases. (4) We 
show results for a 16-point finite-element approximation for dBl./dx^ which 
approximates the derivative at using measurements in a four pixel by four 
pixel region bounded by i — 1 and i + 2 in the x'-direction and j — 1 and j + 2 in 
the y'-direction. Figure 5 shows that the results produced by both this 16-point 
approximation and the three-point approximation are similar. When the global 
minimisation method uses the 16-point finite-element approximation for hori- 
zontal heliographic derivatives we find that the ambiguity-resolution results are 
similar to those retrieved using the three-point approximation, but the 16-point 
approximation requires more iterations for the global minimisation method to 
converge. (5) There are techniques for reliably estimating derivatives from data 
with noise, which produce encouraging results (e.g., Cullum, 1971; Ahnert and 
Abel, 2007). However, there are several challenges in applying these techniques 
to the global minimisation method and we do not attempt to implement these 
techniques in this article. 

5.2. Smoothing the Data Before Disambiguation 

One straightforward approach that may be used to reduce the influence of ran- 
dom noise (like photon noise) is to de-noise or smooth the magnetogram data 
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Figure 5. (a) The absolute value of the difference between the approximations for the 
derivative dBlj./dx^ applied to the no-noise and low- noise ambiguity resolved answer for the 
multipole field configuration, as a function of and j/' (underlying grey-scale); the range 
shown is [10,60] G/pixel. The approximation used for dB]j./dx^ is a three-point finite-dif- 
ference approximation (for details see text). The contours are the same as in Figure 3. (b) 
Same as (a) except the approximation for dB'^/dx^ is a 16-point finite-element approximation 
(for details see text), (c) B^^ as a function of x' at y' = 50 pixels. The black/grey curve 
is for the no-noise/low-noise case; these curves are indistinguishable at many of the points 
shown, (d) The absolute value of the difference between the two curves shown in (c). (e) The 
three-point finite-difference approximation for dB'^/dx^ (for details see text) as a function of 
x' at J/' = 50 pixels. The black/grey curve is for the no-noise/low-noise case, (f) The absolute 
value of the difference between the two curves shown in (e). (g) and (h) Same as (e) and (f) 
except the approximation for dBl,/dx^ is a 16-point finite-element approximation (for details 
see text). 

before disambiguation. To test the efficacy of smoothing the data, we take the 
foUowing approach: (1) Set the initial configuration of azimuthal angles such 
that By > at each pixel. This is done to avoid situations where neighbouring 
pixels get averaged together that initially have very different azimuthal angles. 
(2) Apply a two-dimensional spatial smoothing operator to all three components 
of the magnetic field, B^., By, B\, at both heights; this operation does not change 
the effective pixel size in any direction. (3) Find the minimum of E (Equation (7)) 
using the algorithm described in Section 4.3. (4) Take the realisation of the 
transverse component of the unsmoothed magnetic field that is closest to the 
retrieved solution at each pixel. 
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We have tested several different smoothing operators, including both average 
(i.e., box-car) and triangular smoothing over areas of 3 x 3 pixels and 5x5 pixels. 
We find that the results retrieved with these different smoothing operators tend 
to be quite similar. For the sake of brevity and to highlight the general trends of 
smoothing the data before disambiguation, we show the results for a smoothing 
operator that takes the two-dimensional box-car average over an area of 3 x 3 
pixels. 

Results for the solutions that produce the lowest value of E are provided in 
Figures 6 and 7 and Tables 1 and 2. For the multipole field positioned away 
from disk centre without photon noise this approach does not retrieve the cor- 
rect solution at every pixel over both heights, indicating that the smoothing 
procedure introduces some errors when noise is absent. We have confirmed that 
this is also the case for both of the noise-free test cases examined in Paper I. 
For the cases with photon noise we find that there is some improvement in the 
retrieved solutions in some of the regions with a weak transverse component 
of the magnetic field around the edges of the field of view (especially for the 
high-noise case). However, there are regions where the incorrect solution is 
retrieved for the case with smoothing where the correct solution is retrieved 
for the case where the data is not smoothed before disambiguation. For the tests 
of limited spatial resolution we find that the results retrieved by the method that 
smoothes the data before disambiguation are generally worse than those from 
the method that does not smooth the data (this problem tends to get worse as 
the number of pixels in the smoothing operator gets larger). Given these mixed 
results, we conclude that the smoothing the data before disambiguation does 
not substantially improve the performance of the global minimisation method. 








50 100 150 200 250 300 50 100 150 200 250 300 50 100 150 200 250 300 



Figure 6. Same as the lower panels of Figure 3, except the magnetogram data are smoothed 
before disambiguation with the global minimisation method. 

5.3. The Global Minimisation Method with Additional Constraints 

Several methods that are closely related to the global minimisation method 
incorporate into the minimisation problem an additional smoothness constraint, 
which is based on the current density (Metcalf, 1994; Metcalf et ai, 2006; Leka, 
Barnes, and Crouch, 2009; Leka et al., 2009). In this section we examine a 
couple of approaches that incorporate different additional constraints into the 
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Figure 7. Same as the lower panels of Figure 4, except the magnetogram data are smoothed 
before disambiguation with the global minimisation method. 

global minimisation method. One of the main reasons for using the divergence- 
free property of magnetic fields to resolve the azimuthal ambiguity is to avoid 
unrealistic assumptions. We acknowledge that incorporating a constraint into 
the minimisation problem (in addition to the divergence of the field) may not 
be not entirely consistent with this philosophy. However, in this section we will 
show that the inclusion of an additional constraint can improve the performance 
of the global minimisation method, especially when photon noise is present in 
the data. 

5.3.1. Additional Constraints Involving the Current Density 

Here we consider an approach that assumes that the correct configuration of az- 
imuthal angles over the field of view is the one that corresponds to the minimum 
of 

Ej^ = £^ [|(V-B),,,,fc| + |(V.B),,,,fc+i| +mA,/, {\{J%.,A + 

(8) 

where Aj^ is a dimensionless parameter, /i is the magnetic permeability, and 
{J^)ij,k is the approximation for the vertical heliographic component of the 
current density at pixel {i,j,k). This approach is similar to that taken by the 
minimum energy method MEO (Leka, Barnes, and Crouch, 2009; Leka et ai, 
2009) , where the vertical current density is included in the minimisation problem 
as an additional constraint. In terms of observable quantities (i.e., derivatives 
of the image components of the field with respect to a;'\ y^, and z'), the exact 
expression for is 

dBl dB'l dBl dBl dB'l dB[ 

which does not require the line-of-sight derivative of any of the components 
of the magnetic field. We approximate the horizontal heliographic derivatives 
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using measurements at one height in the same manner as in (V-B)jj-^fc. Conse- 
quently, the approximation ( )i only depends on measurements at the height 
corresponding to the index k. 

The optimal value of parameter Aj, may be problem dependent. We have 
tested several different values of Aj^ in the range from Aj^ = to Aj_ — 10. 
For the set of available test cases we find that values for Aj^ of approximately 
unity consistently produce the best results; we use Aj, = 1 in what follows. To 
find the configuration of azimuthal angles corresponding to the minimum of Ej_^ 
we use the same simulated annealing algorithm that was used for the global 
minimisation method described in Section 4.3. Results for the solutions that 
produce the lowest value of Ej^ are provided in Figures 8 and 9 and Tables 1 
and 2. 




Figure 8. Same as the lower panels of Figure 3 except Ej^ (see Equation (8)) is minimised. 

For the multipole field positioned away from disk centre without photon noise 
this approach retrieves the correct solution at every pixel over both heights; we 
have confirmed that this is also the case for both of the noise-free test cases 
examined in Paper I. As expected, as the level of photon noise increases the 
quality of the results produced by this method decreases (Figure 8 and Table 1). 
However, for the noise-added cases the results retrieved by this approach are 
better than those retrieved by the global minimisation method that minimises 
|V-B| only. This is most pronounced in the regions with a weaker transverse 
component of the magnetic field (for example, note the trend in jM_b^>ioo g in 
Table 1; see also Figure 14). For the noise-added cases we find that Ej^ retrieved 
by the minimisation algorithm is less than that for the answer, indicating that 
the assumption made by this approach is not correct for these test cases over 
the entire field of view. 

For the tests of limited spatial resolution the underlying magnetic field is 
potential, and, therefore, minimising the current density is somewhat justified 
in this case. However, calculation of the vertical component of the current density 
requires the approximation of horizontal heliographic derivatives. Consequently, 
unresolved structure can pose a challenge for this approach (Figure 9 and Ta- 
ble 2). In the plage region at higher spatial resolution we find that the results are 
slightly better than those from the global minimisation method that minimises 
I V-B| only. Away from the plage region, for the higher resolution test cases, there 
is not much difference between the results retrieved by the two approaches. At 
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low resolution (pixel size 0.9") in and around areas with significant unresolved 
structure (see Figure 2) the results retrieved by this approach are better than 
those retrieved for the approach that minimises |V-B| only. For the test cases 
with pixel sizes 0.15", 0.3" and 0.9" we again find that Ej^ retrieved by the 
minimisation algorithm is less than that for the answer. 




100 200 300 400 500 600 



50 100 150 200 250 300 



20 40 60 80 100 



Figure 9. Same as the lower panels of Figure 4 exeept Ej^ (see Equation (8)) is minimised. 

We now consider an approach with an additional constraint that involves all 
three components of the current density vector (Metcalf, 1994; Metcalf et ai, 
2006). In this case, the correct configuration of azimuthal angles over the field 
of view is assumed to be the one that corresponds to the minimum of 



Ej = 5^5^ [|(V-B),j-fe| + |(V-B),,,-fc+i| + Aj/i(J,,,,fc + J.,j,fe+i)] (9) 
i=i j=i 



where Aj is a dimensionless parameter (we use Aj = 1 ), Jij^k = ({Jx)! j ki'^'^y )i j k 

with {J}^)ij^k and {Jy)ij^k the approximations for the x^- and y^-components of 
the current density at pixel (i, j, k), respectively. In terms of observable quantities 
{i.e., derivatives of the image components of the field with respect to x^, and 
z'), the exact expressions for these components of the current density are 
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-aiia23 



(013012 



,9Bl 
032033 j 



+ "33) 



dBi 



013023 



These equations show that the computation of and Jy generally requires the 
line-of-sight derivatives of all three components of the magnetic field, depending 
on the values of an, ai2, 021 and 022- As stated previously, the line-of-sight 
derivatives of the transverse components of the field, and dB'y/dz\ 

depend on the ambiguity resolution at all heights used to approximate line-of- 
sight derivatives. Consequently, for ambiguity-resolution algorithms based on 
the and Jy, the ambiguity resolution at one height affects the ambiguity 
resolution at the other heights used to approximate these derivatives. 

Results for the solutions that produce the lowest value of Ej are provided in 
Figures 10 and 11 and Tables 1 and 2. Broadly speaking, results produced by 
minimising Ej are very similar to those produced by minimising Ej^ . We again 
find that Ej retrieved by the minimisation algorithm is less than that for the 
answer. For the test cases with photon noise the results for the approach that 
minimises Ej are slightly worse than those for Ej^ (see Figure 14). For the tests 
of limited spatial resolution, the results for the approach that minimises Ej are 
slightly better than those for Ej^ for the higher resolution cases, predominantly 
in the plage region. For the lowest resolution case (with pixel size 0.9") the 
results for the approach that minimises Ej are slightly worse than those for 
Ej^ , again, predominantly in the plage region. 




50 100 150 200 250 300 



50 100 150 200 250 300 




Figure 10. Same as Figure 8 except Ej (see Equation (9)) is minimised. 



5.3.2. An Additional Constraint Based on Smoothness of the Magnetic Field 

In this section we consider an approach where the additional constraint included 
in the minimisation problem is designed to suppress pixel-to-pixel variations in 
the magnetic field. In this case the correct configuration of azimuthal angles over 
the field of view is assumed to be the one that corresponds to the minimum of 

= ££(l(V-B),,„fe| + |(V-B),,,-fc+i| + A,^,,,,fc) , (10) 
j=i j=i 
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Figure 11. Same as Figure 9 except Ej (see Equation (9)) is minimised. 

where As is a dimcnsionless parameter and Sij^k is defined as 

S^,,,k = { [4 J, fc) - 4 + 1, k)] ' + [4 (z, J, k) - 4 + 1, A;)] ' }'^' 
+ ^ { [4 (z, J, fc) - i?i J + 1, fc)] ' + [4 (z, j, k) - 4 (z, J + 1, k)]"]"^ 

+ { [4 («, j, fc + 1) - 4 (z + 1, J, fc + 1)] ' + [4 {i, j, fc + 1) - 4 (z + 1, J, fc + 1)] '} 

+ ^ { [4 J, fc + 1) - 4 (z, j + 1, fc + 1)] ' + [4 (z, J, A; + 1) - 4 (z, J + 1, A; + 1)] ' } 

+ ^{[4(*,j',fc)-4(*,j,fc + i)]' + [4(*'^'^)-4(*'j'fc + i)]'}'^' ' (11) 

where di is the distance between the centres of the pixels at and (z + 
at fixed k, and d2 is the distance between the centres of the pixels at and 
(z, J + 1) at fixed k. At the boundaries of the field of view i = Ux and j = rzj, the 
definition of S'i j-,fc is modified to use only pixels within the field of view, in the 
same manner as the finite differencing stencil used to approximate horizontal 
heliographic derivatives. The purpose of S'ij-,fc is to minimise the difference be- 
tween the magnetic field in neighbouring pixels, with the difference normalised 
by the physical distance between measurements. It is worth noting that S'ij.fc 
includes terms that depend on the ambiguity resolution at both of the heights 
k and fc + 1 and, therefore, the ambiguity resolution produced by this approach 
at one height depends on the ambiguity resolution at the other height used in 
the comparison. 

We find that the optimal value of parameter As is problem dependent. We have 
tested several values for As in the range As = to As = 20. We find that values of 
approximately As = 1 consistently produce reasonable results for the various test 
cases (notable exceptions are discussed below); we set As = 1 in what follows. 
The simulated annealing algorithm described in Section 4.3 is used to search for 
the configuration of azimuthal angles that corresponds to the minimum of Eg. 
Results for the solutions that produce the smallest value of found are given 
in Figures 12 and 13 and Tables 1 and 2. 
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For the multipole field configuration without noise this approach retrieves 
the correct solution at each pixel over both heights: This is also the case for 
both of the noise-free test cases examined in Paper I. For both noise-added test 
cases this approach produces better results than both the global minimisation 
method that minimises |V-B| only and that which minimises a combination of 
|V-B| and the current density {i.e., Ej^ or Ej). This is especially noticeable in 
regions with x' w 0, a;' w 300, and at w 170 and y' « 50, where the transverse 
component of the magnetic field is relatively weak and the signal tends to be 
dominated by noise (see Figure 1 and Figure 14). Again, for the noise-added test 
cases we find that E^ retrieved by the minimisation algorithm is less than that 
for the answer. For the noise-added test cases we also find that slightly better 
results than those shown in Figure 12 can be achieved with this approach in the 
regions with a weaker transverse component of the field with values of Ag slightly 
larger than unity. 

Broadly speaking, for the tests of limited spatial resolution the results pro- 
duced by this approach are not much different to those retrieved by the global 
minimisation method that minimises |V-B| only and that which minimises a 
combination of | V-B| and the current density. There are some subtle differences 
that are worth noting for the case with pixel size 0.9". The results produced by 
minimising Eg tend to be a little better in the plage region than those produced 
by minimising E or Ej^, but a little worse around the flux concentration at 
the top, right of the field of view. Nevertheless, it is encouraging that this 
approach performs relatively well for the tests of limited spatial resolution, as the 
assumption that the magnetic field is smooth made by this approach may not be 
correct in regions with significant unresolved structure. Again, for the test cases 
with pixel sizes 0.15", 0.3" and 0.9" we find that E^ retrieved by the minimisation 
algorithm is less than that for the answer, indicating that the assumption made 
by this approach is indeed not correct. For the tests of limited spatial resolution 
we find that the results produced by minimising Eg tend to get slightly worse 
for values of Ag larger than unity, especially for the case with pixel size 0.9". 

We have also tried an approach where the additional smoothness constraint 
in the minimisation problem is based on the second horizontal heliographic 
derivatives of the magnetic field (results not shown). For the tests of Poisson 
noise the results produced by this approach are similar to those produced by 
minimising Ej^ but not as good as the results produced by minimising Eg. For 
the tests of limited spatial resolution the results produced by the approach based 
on second derivatives are similar to those produced by minimising Ej^ or Eg. 

5.4. A Hybrid Method: the Global Minimisation Method Combined with a 
Smoothing Algorithm 

Evidently, the global minimisation method that seeks to minimise a combination 
of the divergence and the diff'erence between the magnetic field in neighbouring 
pixels can produce better results than the global minimisation method based only 
on the divergence. However, none of the various global minimisation methods im- 
plemented in this article produce desirable results in regions where the transverse 
component of the magnetic field is strongly aff'ected by photon noise. Figure 14 
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Figure 12. Same as Figure 8 except Eb (see Equation (10)) is minimised. 




Figure 13. Same as Figure 9 except Eb (see Equation (10)) is minimised. 



shows the fraction of pixels correctly resolved as a function of the magnitude of 
the transverse component of the magnetic field for the two test cases with photon 
noise. This figure clarifies how the various methods perform in regions with a 
weaker transverse component of the magnetic field where the influence of photon 
noise is most pronounced. Figure 14 shows that for each method as B± decreases 
the fraction of pixels correctly resolved generally decreases, for low values of the 
transverse component of the magnetic field. Figure 14 also shows that, of the 
global minimisation methods, the one that minimises | V-B| only performs worst 
and the one that minimises Eg performs best. The global minimisation methods 
with an additional constraint involving the current density produce intermediate 
results for these test cases. We note that in practice the choice for the additional 
constraint used by the global minimisation method may need to be adjusted 
depending on the expected nature of the field. 

Given the results for the global minimisation methods shown in Figure 14 we 
present a hybrid method that first minimises Eg (Equation (10)), as described 
in Section 5.3.2, then revisits pixels in regions where the transverse component 
of the magnetic field is weak (determined by a prescribed "threshold") with 
a smoothing algorithm that does not use the divergence-free condition. The 
smoothing algorithm implemented here is similar to that used by MEO (Leka, 
Barnes, and Crouch, 2009; Leka et ai, 2009) and aims to minimise the difference 
between the magnetic field at a target pixel with that in a small neighbourhood 
of surrounding pixels. 
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Figure 14. (a) / (Bj.), the fraction of pixels correctly resolved in the interval [B^, B± + 5B±) 
(where SB± = 20G) as a function of the magnitude of the transverse component of the 
magnetic field, B±, for the test case with a low level of photon noise. The black circles joined 
by a dashed line are for the global minimisation method that minimises |V-B| only. The grey 
circles joined by a full line are for the global minimisation method that minimises Ej^ (see 
Equation (8)). The grey squares joined by a dashed line axe for the global minimisation method 
that minimises Ej (see Equation (9)). The black squares joined by a dashed line are for the 
global minimisation method that minimises Es (see Equation (10)). The black circles joined by 
a full line are for the hybrid method (see Section 5.4). (b) Same as (a) except for the high-noise 
test case. Note that, in both panels, the curves joining the points are only included as a guide. 

The smoothing algorithm proceeds as follows. The azimuthal angles in pixels 
considered above the prescribed threshold are fixed. Pixels below threshold are 
visited according to the following criteria: (1) Pixels with more neighbours con- 
sidered above threshold are visited first. (2) For pixels with the same number of 
neighbours considered above threshold, those with a stronger transverse compo- 
nent of the magnetic field are visited first. Once a pixel has been examined it is 
treated as if it is above threshold. After each pixel is visited the neighbourhood 
information for each below-threshold pixel surrounding it is updated to ensure 
that pixels arc visited in the order described above. At each target pixel we select 
the configuration of azimuthal angles that corresponds to the smallest average 
difi^erence between the magnetic field in the target pixel and that in surround- 
ing pixels. In doing so, the azimuthal angles in any of the surrounding pixels 
considered below threshold are allowed to change independently. Consequently, 
the results produced by this smoothing algorithm do not depend on the initial 
configuration of azimiithal angles for below-threshold pixels. 

We have experimented with several different definitions for the neighbourhood 
used by the smoothing algorithm. Broadly speaking, we find that neighbour- 
hoods involving more pixels perform better than those with fewer pixels, and 
those that involve pixels from two heights perform better than those that use 
only pixels at one height. In Figures 15 and 16 and Tables 1 and 2 we show results 
for this hybrid approach where the smoothing algorithm uses a neighbourhood 
of five pixels by five pixels at the same height as the target pixel plus the single 
pixel directly above (or below) the target pixel. 

Wc find that the ambiguity-resolution results depend on the prescribed thresh- 
old. For the test cases with photon noise, the ambiguity-resolution results gen- 
erally improve as the threshold increases from zero to a point whereafter the 
results gradually get worse. This indicates that in practice the threshold may 
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need to be adjusted according to the level of photon noise present in the data. 
For demonstration purposes we use a threshold of 150 G for the low-noise case 
and 400 G for the high-noise case: These thresholds arc chosen because they cor- 
respond roughly to the value where the smoothing algorithm produces optimal 
results in experiments. In these cases, the smoothing algorithm clearly produces 
results that are greatly improved in the regions below the prescribed threshold, 
in comparison to the results produced by minimising Eg (see Figure 14, Figure 15 
and Table 1). 

For the sake of comparison, for the test cases without photon noise we use a 
threshold of 100 G. This value is somewhat arbitrary since there is no photon 
noise in these cases. For cases without photon noise we find that the results 
produced by the hybrid approach are very similar to those produced by min- 
imising Es- However, a small number of pixels below threshold do have their 
azimuthal angles altered by the smoothing algorithm. This indicates that the 
smoothing algorithm may not be necessary, and may be counterproductive, for 
regions where the transverse component of the magnetic is only weakly affected 
by photon noise. 




50 100 150 200 250 300 50 100 150 200 250 300 50 100 150 200 250 300 



Figure 15. Same as Figure 8 except tlie liybrid method is employed, as described in Section 5.4 
(i.e., , first Eb (see Equation (10)) is minimised, then a smoothing algorithm is applied to pixels 
below a threshold transverse magnetic field strength). 




Figure 16. Same as Figure 9 except the hybrid method is employed, as described in Section 5.4 
{i.e., , first Eb (see Equation (10)) is minimised, then a smoothing algorithm is applied to pixels 
below a threshold transverse magnetic field strength). 
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Figure 14 shows that a reasonable hybrid method may also be constructed 
by combining one of the other global minimisation methods with the smoothing 
algorithm, although this may require a higher threshold (ideally, the threshold 
for the smoothing algorithm is chosen such that all pixels above the threshold 
are correctly resolved). For example, in other experiments we apply the smooth- 
ing algorithm to the results produced by the global minimisation method that 
minimises |V-B| only. In this instance, to obtain optimal solutions for the low- 
noisc test case we find that the threshold needs to bo increased to approximately 
250 G; the results produced are comparable to those shown in Figure 15. We 
acknowledge that it may not be desirable in practice to use the smoothing 
algorithm over large areas in general because it does not use the divergence- 
free condition. If this is the case, the results in Figure 14 suggest that the most 
reliable hybrid method is the one that first uses the global minimisation method 
that minimises Eg, as this allows the smoothing algorithm to start from a lower 
threshold. 

6. Conclusions 

Wc investigate how the azimuthal ambiguity, present in solar vector magne- 
togram data, can be resolved by using the divergence-free property of magnetic 
fields, along with information regarding the variation of the field in the line-of- 
sight and horizontal heliographic directions. In Paper I several methods were 
tested that each make different assumptions about how to use the divergence 
to resolve the ambiguity, using error-free synthetic data that did not account 
for noise that is typical of solar observational data. In this article we use more 
realistic synthetic data to test how the various methods examined in Paper I 
respond to the effects of two kinds of noise: noise to simulate Poisson photon 
noise in the observed polarization spectra, and a spatial binning to simulate the 
effects of limited instrumental spatial resolution in the directions perpendicular 
to the line-of-sight. We find the same general trend that was found in Paper I 
in that the global minimisation method is more robust than both the Wu and 
Ai criterion and the sequential minimisation method (Boulmezaoud and Amari, 
1999). However, we find that all methods based on the divergence-free property 
examined in Paper I can produce undesirable results when photon noise or 
unresolved structure are present in the data. 

We conduct a series of experiments aimed at improving the performance of the 
global minimisation method and reducing the effects of noise on the disambigua- 
tion results. Based on the results from these experiments we present a two-step, 
hybrid method that produces reasonable results in tests with synthetic data. The 
first step of the hybrid method is similar to the global minimisation method of 
Paper I and involves the minimisation of a combination of the approximation for 
the divergence and a smoothness constraint that seeks to minimise the difference 
between the magnetic field in neighbouring pixels. We test several versions of this 
modified global minimisation method that have different additional constraints 
and find that all of the methods tested produce undesirable results in regions 
where the transverse component of the magnetic field is strongly affected by 
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photon noise. Consequently, in the second step of the hybrid approach pixels 
with measurements of the transverse component of the magnetic field below a 
prescribed threshold are revisited with a smoothing algorithm that also seeks to 
minimise the difference between the magnetic field in neighbouring pixels. For 
synthetic data with photon noise the hybrid approach produces results that are 
significantly better than those produced by the global minimisation method that 
minimises only the divergence. For synthetic data with unresolved structure the 
hybrid approach produces results that are slightly bcittcr than those produced by 
the global minimisation method that minimises only the divergence. The hybrid 
method uses an assumption about the smoothness of the magnetic field that is 
not explicitly related to the divergence of the field. We acknowledge that this 
assumption may not be appropriate for solar magnetic fields in general. 

For disambiguation methods based on the divergence-free condition there are 
some sources of uncertainty that are not modelled by the synthetic data employed 
in this investigation. First, limited spatial resolution in the line-of-sight direction; 
the synthetic data employed here assume that the locations of the observation 
heights arc known exactly, but this may not be accurate for solar observational 
data. Second, the synthetic data is constructed assuming that the magnetic field 
is measured at constant geometrical height over the field of view, but in practice 
this may not be a reasonable approximation as the solar atmosphere may vary 
significantly in the directions perpendicular to the line-of-sight. Third, the line- 
of-sight distance between the two observation heights is assumed to be constant, 
which may not be a reasonable approximation for the similar reasons. These 
additional sources of uncertainty may be significant for solar observational data 
and further investigation is required to quantify how these issues may affect the 
results produced by the methods tested in this article. 

In the interest of further development, testing and implementation, the soft- 
ware for the global minimisation methods and the hybrid method is available at 
http: / /www.cora. nwra .com /'^ash / ambig2.tar.gz. 

Acknowledgements The author thanks K.D. Leka and Graham Barnes for providing the 
synthetic datasets used in this investigation. Tlic author thanks the referee for lielpful sugges- 
tions. This work was supported by funding from NASA under contracts NNH05CC49C/NNH05CC75C| 
and NNH09CE60C. This work used the Extreme Science and Engineering Discovery Envi- 
ronment (XSEDE), which is supported by National Science Foundation grant number OCI- 
1053575. 



References 

Ahnert, K., Abel, M.: 2007, Comput. Phys. Commun. 177, 764. doi:10.1016/j.cpc.2007.03.009. 

Boulmezaoud, T.Z., Amari, T.: 1999, Astron. Astrophys. 347, 1005. 

Chiu, Y.T., Hilton, H.H.: 1977, Astrophys. J. 212, 873. doi:10.1086/155111. 

CoUados, M., Martinez Fillet, V., Ruiz Cobo, B., del Toro Iniesta, J.C., Vazquez, M.: 1994, 

Astron. Astrophys. 291, 622. 
Crouch, A.D., Barnes, G.: 2008, Solar Phys. 247, 25. doi:10.1007/sll207-007-9096-l. 
Grouch, A.D., Barnes, G., Leka, K.D.: 2009, Solar Phys. 260, 271 (Faper I). 

doi:10.1007/sll207-009-9454-2. 
GuUum, J.: 1971, SIAM J. Numeric. Anal. 8, 254. doi:10.1137/0708026. 



ms.tex; 11 October 2012; 0:23; p. 28 



Rcsohiug llw Aziiimilial Aml)iguit\- iu Vector Maguotogram Data 



29 



Cuperman, S., Li, J., Semel, M.: 1993, Astron. Astrophys. 278, 279. 

del Toro Inicsta, J.C., Ruiz Cobo, B.; 1996, Solar Phys. 164, 169. doi:10.1007/BF00145631. 
Eibc, M.T., Aulanier, G., Faurobcrt, M., Mcin, P., Maliicrbc, J.M.: 2002, Astron. Astrophys. 

381, 290. doi:10.1051/0004-5361:20011495. 
Gary, G.A., Hagyard, M.J.: 1990, Solar Phys. 126, 21. doi:10.1007/BF00158295. 
Gcorgoulis, M.K.: 2012, Solar Phys. 276, 423. doi;10.1007/sll207-011-9819-l. 
Harvey, ,I.W.: 1969, Magnetic Fields Associated with Solar Active-Region Prominences. Ph.D. 

thesis, University of Colorado, Boulder. 
Kirkpatrick, S., Gelatt, C.D., Vecchi, M.P.: 1983, Science 220, 671. 

doi:10.1126/science.220. 4598.671. 
Leka, K.D., Barnes, G.: 2012, Solar Phys. 277, 89. doi;10.1007/sll207-011-9821-7. 
Leka, K.D., Metcalf, T.R.: 2003, Solar Phys. 212, 361. doi:10.1023/A:1022995404064. 
Leka, K.D., Barnes, G., Crouch, A.: 2009, In; Lites, B., Cheung, M., Magara, T., Mariska, 

J., Reeves, K. (eds.) The Second Hinode Science Meeting: Beyond Discovery-Toward 

Understanding, ASP Conf. Ser. 415, 365. 
Leka, K.D., Barnes, G., Crouch, A.D., Metcalf, T.R., Gary, G.A., Jing, J., Liu, Y.: 2009, Solar 

Phys. 260, 83. doi:10.1007/sll207-009-9440-8. 
Leka, K.D., Barnes, G., Gary, G.A., Crouch, A.D., Liu, Y.: 2012, Solar Phys. 276, 441. 

doi;10.1007/sll207-011-9879-2. 
Li, J., Amari, T., Fan, Y.: 2007, Astrophys. J. 654, 675. doiilO. 1086/509062. 
Li, J., Cuperman, S., Semel, M.: 1993, Astron. Astrophys. 279, 214. 
Liu, Y., Wang, J., Yan, Y., Ai, G.: 1996, Solar Phys. 169, 79. doi:10.1007/BF00153835. 
Maltby, P., Avrett, E.H., Carlsson, M., Kjeldseth-Moe, O., Kurucz, R.L., Loeser, R.: 1986, 

Astrophys. .J. 306, 284. doiilO. 1085/154342. 
Metcalf, T.R.: 1994, Solar Phys. 155, 235. doi:10.1007/BF00680593. 

Metcalf, T.R., Jiao, L., McClymont, A.N., Canfield, R.C., Uitenbroek, H.: 1995, Astrophys. 

J. 439, 474. doi:10.1086/175188. 
Metcalf, T.R., Leka, K.D., Barnes, G., Lites, B.W., Georgoulis, M.K., Pevtsov, A. A., Bala- 

subramaniam, K.S., Gary, G.A., Jing, J., Li, J., Liu, Y., Wang, H.N., Abramenko, V., 

Yurchyshyn, V., Moon, Y.J.: 2006, Solar Phys. 237, 267. doi:10.1007/sll207-006-0170-x. 
Metropolis, N., Rosenbluth, A.W., Rosenbluth, M.N., Teller, A.H., Teller, E.: 1953, J. Chem. 

Phys. 21, 1087. doi:10.1063/l. 1699114. 
Press, W.H., Teukolsky, S.A., Vetterling, W.T., Flannery, B.P.: 1992, Numerical Recipes in 

FORTRAN. The Art of Scientific Computing, Cambridge University Press, Cambridge, 

436. 

Ruiz Cobo, B., del Toro Iniesta, J.C.: 1992, Astrophys. J. 398, 375. doi: 10. 1086/171862. 
Socas-Navarro, H.: 2005, Astrophys. J. Lett. 631, L167. doi:10.1086/497334. 
Socas-Navarro, H.: 2007, Astrophys. J. Suppl. 169, 439. doi:10.1086/510336. 
Socas-Navarro, H., Trujillo Bueno, J., Ruiz Cobo, B.: 2000, Astrophys. J. 530, 977. 
doi:10.1086/308414. 

Vernazza, J.E., Avrett, E.H., Loeser, R.: 1981, Astrophys. J. Suppl. 45, 635. 
doi:10.1086/190731. 

Westendorp Plaza, C, del Toro Iniesta, ,I.C., Ruiz Cobo, B., Martinez Pillet, V., Lites, B.W., 
Skumanich, A.: 1998, Astrophys. J. 494, 453. doiilO. 1085/305192. 

Westendorp Plaza, C, del Toro Iniesta, ,I.C., Ruiz Cobo, B., Martinez Pillet, V., Lites, B.W., 
Skumanich, A.: 2001, Astrophys. J. 547, 1130. doi:10.1086/318376. 

Wu, L.X., Ai, G.X.: 1990, Acta Astrophys. Sinica 10, 371. 



ms.tex; 11 October 2012; 0:23; p. 29 



ms.tex; 11 October 2012; 0:23; p. 30 



